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(54) High resolution 3D ultrasound imaging system deploying a multi-dimensional array of 
sensors and method for multi-dimensional beamforming sensor signals 



(57) An adaptive multidimensional beamformer 
liaving near-instantaneous convergence for ultrasound 
imaging systems deploying multidimensional sensor ar- 
rays is disclosed. In a first step, the multidimensional 
beamfomier is decomposed into sub-apertures. Each 
sub-aperture is then again decomposed into two coher- 
ent subsets of circular and/or line array beamfomiers in 
different coordinate directions of the multidimensional 
array. Implementation of the multidimensional beam- 



former according to the present invention provides the 
basis for a 3D ultrasound imaging system according to 
the present invention comprising a compact multidimen- 
sional sensor array and a compact processing unit that 
Is field deployable and generates high resolution three- 
dimensional images in real time. It is also possible to 
capture four-dimensional images, the fourth dimension 
being time and the resulting images fomning a video im- 
age of a volume of a moving organ. 
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Description 

Field of the Invention 

5 [0001] This Invention relates to the field of ultrasound imaging systems. In particular, the present invention relates 
to field deployable 3D ultrasound Imaging systems providing high resolution images In real time. 

Baclcground of the Invention 

10 [0002] Ultrasound sensing and imaging technology provides a powerful tool for non-Invasive Imaging for treatment 
assessment and for minimally Invasive surgery. Unlike CT scanners and MRI ultrasound Imaging systems are compact 
and much cheaper to manufacture. These advantages allow use of ultrasound imaging systems in mobile units such 
as an ambulance or a helicopter. In general, victims of accidents, disasters or wars need immediate assessment and 
treatment in order to save their lives. For example, deployment of compact ultrasound imaging systems In mobile units 

15 allows on site imaging for treatment assessment during transportation providing live saving information for later surgery 
in a hospital or even providing information for minimally Invasive surgery within the mobile unit. Therefore, It would be 
highly advantageous to provide a compact field-depioyable 3D ultrasound Imaging system for mobile units and field 
hospitals for immediate Imaging of victims on accident sites, in disaster areas or in war zones. 
[0003] However, state of the art ultrasound Imaging systems suffer from very poor image resolution due to the very 

20 small size of sensor arrays of compact systems. Therefore, such systems do not provide Images having a satisfying 
resolution for treatment assessment or surgery. In order to improve image quality It Is necessary to deploy a large 
number of sensors In a compact multidimensional array to provide significant Improvements in array gain for signals 
embedded in partially correlated noise fields. Partially correlated noise fields are caused, for example, by non-linear 
propagation characteristics of the human body and result in aberration effects and fuzziness In reconstructed Images. 

25 The improvements in array gain result in image resolution improvements and minimization of the aberration effects. 
[0004] An overview of the state of the art in adaptive and synthetic aperture beamfonners is given in " Implementation 
of Adaptive and Synthetic Aperture Processing Schemes in Integrated Active- Passive Sonar Systems Proceedings 
of the IEEE, 86(2), pp. 358-397, February, 1998 by S. Stergiopoulos. These algorithms have been designed to increase 
the signal - to • noise ratio for improved target detection and to provide simultaneously parameter estimates such as 

30 frequency, time delay, Doppler shift and bearing for Incorporation Into algorithms localising, classifying and tracking 
acoustic signals. 

[0005] To optimize the beam forming process, beamfomning filter coefficients have to be chosen based on data 
received from a sensor array of the sonar system. I n particular, the coefficients have to be chosen based on the statistical 
properties of the received data. Algorithms using characteristics of noise received from the sensor array for optimizing 

35 the beamforming process are called adaptive beamformers. The adaptive beamformers require knowledge of a cor- 
related noise's covariance matrix. However, If the knowledge of the noise's characteristic Is Inaccurate, performance 
of the adaptive beamformer will degrade significantly and may even result in cancellation of a desired signal. Therefore, 
it Is very difficult to Implement useful adaptive beamformers in real time operational systems. Numerous articles on 
adaptive beamfomners and the difficulties concerning their implementation have been published. Various adaptive 

40 beamformers such as the Generalized SIdelobe Cancellers (GSC), the Linearly Constrained Minimum Variance Beam- 
formers (LCMV) and the Minimum Variance Distortionless Response (MVDR) are discussed In the following references, 
which are hereby incorporated by reference: 

B. Windrow et al.: "Adaptive Antenna Systems", Proceedings IEEE, 55(12), pp. 2143-2159, 1967; 
45 N.L. Owsley: "Sonar Array Processing", S. Haykin, Editor, Prentice-Hall Signal Processing Series, A.V. Oppenheim 

Series Editor, pp. 123, 1985; 

B. Van Veen and K. Buckley: "Beamforming: a Versatile Approach to Spatial Filtering", IEEE ASSP Mag., pp. 4-24, 
1988; 

J. Capon: "High Resolution Frequency Wavenumber Spectral Analysis", Proc. IEEE, 57, pp. 1408-1418, 1969; 
50 s. Haykin: "Adaptive Filter Theory", Prentice-Hall, Englewood Cliffs, NJ, 1986; 

TL. Marzetta: "A New Interpretation for Capon's Maximum Likelihood Method of Frequency-Waven umber Spectra 

Estimation", IEEE-Trans. Acoustic Speech Signal Proc, ASSP-31(2), pp. 445-449, 1983; 

A.H. Sayed and T. Kailath: "A State-Space Approach to Adaptive RLS Filtering", IEEE SP Mag., pp. 18-60, July, 

1994; 

55 A.B. Baggeroer, W.A. Kuperman and PN, Mikhalevsky: "An Overview of Matched Field Methods in Ocean Acous- 

tics", IEEE J. Oceanic Eng., 18(4). pp. 401-424. 1993; 

H, Wang and M. Kaveh: "Coherent Signal-Subspace Processing for the Detection and Estimation of angles of 
Arrival of Multiple Wideband Sources", IEEE Trans. Acoust. Speech, Signal Proc, ASSP-33, pp. 823-831 , 1985; 
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J. Krolik and D.N. Swingler: "Bearing Estimation of Multiple Brodband Sources using Steered Covariance jVlatrices", 
IEEE Trans. Acoust. Speech, Signal Proc, ASSP-37, pp. 1481-1494, 1989; 

S.D. Peters: "Near-Instantaneous Convergence for Memoryless Narrowband GSC/NLMS Adaptive Beamformers", 
submitted to IEEE Trans. Acoust. Speech, Signal Proc, January 1995; 
5 L J. Griffiths and C.W. Jim: "An Alternative Approach to Linearly Constrained Adaptive Beamforming", IEEE Trans, 

on Antennas and Propagation, AP-30, pp. 27-34, 1982; and, 

D.T.M. Slock: "On the Convergence Behavior of the LMS and the Nomnalized LMS Algorithms", IEEE Trans. Acoust. 
Speech, Signal Proc, ASSP-31, pp. 2811-2825, 1993. 

10 [0006] Unfortunately, implementation of adaptive beamformers in modern ultrasound systems comprising multi-di- 
mensional arrays with hundreds of sensors requires very large amounts of memory and very large processing capa- 
bilities for real time data processing making their application for field-deployable systems impossible. To implement 
adaptive beamfomners using current computer technology, the concept of partially adaptive beamfonner design has 
been developed. The partially adaptive beamfonner reduces the number of degrees of freedom, associated with the 

15 beamforming process, lowering the computational requirements and improving response time. Unfortunately, due to 
the reduction of the number of degrees of freedom the partially adaptive beamformers cannot converge to an optimum 
solution as fully adaptive beamformers. Therefore, the partially adaptive beamformers cannot make substantial use of 
multidimensional an-ays deployed in ultrasound systems in order to improve array gain and to provide images with high 
resolution. 

20 [0007] It is, therefore, an object of the invention to overcome the problems associated with the implementation of 
adaptive beamfonners in modem ultrasound imaging systems comprising multidimensional sensor arrays. 
[0008] It is further an object of the invention to provide adaptive beamformers with near-instantaneous convergence 
for ultrasound imaging systems deploying line arrays, circular arrays, spherical arrays of sensors or any superposition 
of these types of arrays. 

25 [0009] It is yet another object of the invention to provide a 3D ultrasound imaging system comprising a multidimen- 
sional sensor array for generating high resolution Images in real time using an adaptive beamfomning process that is 
field-deployable. 

Summary of the Invention 

30 

[0010] In accordance with the present invention there is provided, an adaptive multidimensional beamformer having 
near-instantaneous convergence for ultrasound imaging systems. Implementation of the multidimensional beamformer 
according to the present invention provides the basis for a 3D ultrasound imaging system according to the present 
invention comprising a compact multidimensional sensor array and a compact processing unit that is field deployable 
35 and generates high resolution images in real time or near real time, 

[001 1 ] In accordance with the present invention there is provided, a method for multidimensional beamforming sensor 
time series provided by sensors deployed in a multidimensional array of an ultrasound imaging system comprising the 
steps of: 

40 decomposing the multidimensional beamformer into two coherent subsets of circular array beamfonners and/or 

line array beamformers, a first subset comprising subsequent beamformers in a first coordinate direction of the 
multidimensional array and a second subset comprising subsequent beamformers in a second other coordinate 
direction of the multidimensional array; 

beamforming for a predetermined beam steering direction of data relating to the sensor time series by applying 
45 the subsequent beamformers of the first subset, each beamfonner producing a beam time series; and, 

beamforming for the predetemnined beam steering direction each beam time series of the first subset of beam- 
formers applying the subsequent beamfomners of the second subset for the steered direction producing one beam 
time series for the beam steering direction, 

so [0012] In accordance with aspect of the present invention there is provided, a method for multidimensional beam- 
forming sensor time series provided by sensors deployed in a multidimensional array of an ultrasound Imaging system 
using a sub-aperture configuration, the method comprising the steps of: 

a) dividing the multidimensional beamformer into a plurality of subsequent sub-apertures; b) decomposing each 
55 sub-aperture into two coherent subsets of circular array beamfonners and/or line array beamfonners, a first subset 

comprising subsequent beamformers in a first coordinate direction of the multidimensional array and a second 
subset comprising subsequent beamformers in a second other coordinate direction of the sub-aperture; 
c) conventional beamforming each sub-aperture for a predetennined beam steering direction a Fourier transform 
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of the sensor time series by applying the subsequent beamformers of the first subset each beamfonner producing 
a beam time series; 

d) conventional beamfomning each sub-aperture for the predetermined beam steering direction the beam time 
series produced by the beamformers of step c) by applying the subsequent beamformers of the second subset for 

5 the steered direction producing one beam time series for the beam steering direction for each sub-aperture; and, 

e) adaptive beamforming on line arrays, each line array comprising beam time series of different subsequent sub- 
apertures in one coordinate direction, providing one or more beam time series for the beam steering direction. 

[0013] In accordance with another aspect of the present invention there is provided, a method for beamfomiing 
10 sensor time series provided by sensors of an ultrasound imaging system using a coherent broad band adaptive beam- 
former, the method comprising the steps of: 

segmenting the continuous sensor time series Into a set of overlapped data sets; 

calculating a FFT of each overlapped data set producing a set of Fourier transforms of the overlapped data sets 
15 for different frequency bins; 

forming a cross spectral density matrix from the Fourier transfomris of the overlapped data sets for each frequency 
bin and each predetennined steering direction; 

fomning a steering covarlance matrix using the cross spectral density matrix and a diagonal matrix of conventional 
steering vectors, one steering covarlance matrix for each steering direction and a frequency band of interest; 
20 inverting the steering covarlance matrices; 

estimating adaptive steering vectors by assuming stationarity across frequency bins of a frequency band of interest 
and considering an estimate of the steering covarlance matrix being the same as a narrow band estimate for a 
center frequency of the frequency band of Interest; 

determining narrow band adaptive steering weights using the estimate of the adaptive steering vectors; 
25 forming adaptive beams in frequency domain from the Fourier transform of the overlapped data sets and the 

adaptive steering weights; 

forming adaptive beams in time domain through IFFT; and, 

detemiining continuous beam time series by discarding overlap and concatenation of segments. 

30 [0014] in accordance the present invention there is further provided, a field deployable 3D ultrasound Imaging system 
for producing high resolution 3D images of an object In real time, the ultrasound system comprising: 

a source for emitting ultrasound waves; 

a compact adaptive multidimensional sensorarray for capturing reflections of the ultrasound waves, the ultrasound 
35 waves being reflected by different structures within the object, and for providing sensor time series indicative of 

the reflected ultrasound waves; 
a compact processing unit for: 

receiving the sensor time series produced by the multidimensional sensor array; 
40 processing the sensor time series in order to produce continuous beamtime series by: 

decomposing a multidimensional beamfonner Into sub-apertures comprising coherent subsets of circular 
array beamformers and/or line array beamfonners; 

conventional beamforming circular arrays; 
45 adaptive beamforming line arrays; 

and, 

reconstructing 3D Images from the beam time series in real time; 

50 and, 

a display for displaying the reconstructed 3D images in real time. 

Brief Description of the Drawings 

55 [001 5] Exemplary embodiments of the invention will now be described in conjunction with the following drawings, in 
which: 

Fig. 1 is a simplified diagram of a compact field deployable ultrasound imaging system according to the invention; 
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Fig. 2 is a simplified diagram of a line array of sensors for an ultrasound imaging system; 
Fig. 3 is a simplified diagram of a circular array of sensors for an ultrasound Imaging system; 
Fig. 4 Is a simplified diagram of a cylindrical array of sensors for an ultrasound Imaging system illustrating a de- 
composition concept according to the invention; 
5 Fig. 5 is a simplified flow diagram of an adaptive beamfomier according to the invention for an ultrasound imaging 

system; 

Fig. 6 Is a simplified diagram of a line array of sensors for an ultrasound imaging system Illustrating an adaptive 
sub-aperture structure according to the Invention; 

Fig. 7 is a simplified diagram of a circular array of sensors for an ultrasound imaging system Illustrating an adaptive 
10 sub-aperture structure according to the invention; 

Fig. 8 is a simplified diagram of a cylindrical array of sensors for an ultrasound imaging system illustrating an 
adaptive sub-aperture structure according to the invention; 

Fig. 9 is a simplified diagram of a planar array of sensors for an ultrasound imaging system illustrating an adaptive 

sub-aperture structure according to the invention; 
IS Fig. 1 0 is a schematic diagram of a signal processing flow according to the invention for beamforming cylindrical 

arrays; 

Fig. 11 Is a schematic diagram of a signal processing for an adaptive sub-aperture structure for multidimensional 
arrays according to the Invention; 

Fig. 12 is a schematic diagram illustrating power of beam response to active CW pulses as function of steering 
20 angle for various types of beamformers; 

Fig. 13 illustrates simulated matched filter results of various types of beamformers for ultrasound applications; 
Fig. 14 is a simplified diagram of an ultrasound Imaging system to assess image resolution improvements of an 
adaptive beamformer according to the invention; and, 

Fig. 15 shows various 3D images of a fetus' skull obtained by using state-of-the-art signal processing techniques 
25 and by applying the adaptive beamformer according to the present invention to the signals captured by a state- 

of-the-art ultrasound system. 

Detailed Description 

30 [0016] Fig. 1 illustrates schematically a field deployable 3D ultrasound imaging system 100 according to the invention. 
An object 2 is irradiated by a source 4 emitting ultrasound waves. Different structures 3 within the object 2 reflect the 
ultrasound waves in different directions. The reflected ultrasound waves are then captured by a plurality of sensors 8 
of a compact multidimensional - 2D or 3D - sensor array 6. Output signals of the sensors 8 are transmitted via a 
communication link to a compact processing unit 12 such as a conventional computer workstation for processing. 

35 Reconstructed images are then displayed using display 14. Provision of high resolution images and compact design 
of the multidimensional sensor array 6 and the processing unit 12 for a field deployable ultrasound imaging system is 
made possible using an advanced beamfonning process according to the present invention which will be disclosed in 
the following sections. The advanced beamforming process according to the present Invention maximizes array gain 
and image resolution by effectively using all information provided by the sensor signals of the multidimensional sensor 

40 array 6. This effective use of the sensor signals allows to minimize the multidimensional sensor array 6, which is 
essential for the field deployable ultrasound imaging system 100. Furthermore, the advanced beamforming process 
according to the invention comprises a decomposition dividing a fully multidimensional beamfomier into sub-sets of 
coherent processes which may be implemented on small size CPU's integrated within a parallel configuration of existing 
computing architectures. Therefore, the advanced beamfonning process according to the present invention is per- 

45 formed using the compact processing unit 12 of the field deployable ultrasound imaging system 100. Additional, use 
of a sub-aperture structure of the beamfomning process provides near-instantaneous convergence of adaptive beam- 
formers allowing real time applications of the ultrasound system according to the present invention. 
[0017] As a result of the decomposition according to the invention the fully multidimensional beamfomning process 
is reduced to beamforming processes for line and circular arrays. 

50 [0018] Fig. 2 illustrates a linear sensor array 20 comprising N sensors 8 with uniform sensor spacing 5 receiving a 
plane wave signal 22 with a direction of propagation 0. 

^«(^iX^ = = 1,...,/) 

55 

are sensor time series produced by the N sensors 8 indicative of the received plane wave signal 22. For the plane 
wave signal 22 arriving at the angle 6 the sensor time series are expressed by 
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15 



20 



5 wherein t(G) is a time delay between the sensor time series of two adjacent sensors given by t(G) = 5cos(G)/c with c 
being a propagation speed of the plane wave signal 22 in a medium of interest. 
[0019] The plane wave signal response of the N - sensor line array 20 is expressed by 

'0 B{/,4) = D'{f,4)x{f), (1) 

wherein * denotes a complex conjugate transpose. X(/) Is a vector of the Fourier transform of 

for a frequency bin / of Interest. 

25 Is a steering vector having, for example, Its n^^ phase temn for the plane wave signal 22 at an angle 0^ being expressed by 

exp[y2^(A7~ l)*os ^/c], 
wherein 6^ is a beam steering direction angle. The beam power pattern 

is then obtained by 

Side lobe structures resulting from the Fourier transformation of discontinuities in the signals, also called Gibbs effect, 
may be suppressed by applying different weights during the Fourier transfomriatlon - called shading - at the expense 
^5 of an increased beam width. Angular response of the line array is ambiguous with respect to the angle 85, responding 
equally to targets at angle 0^ and - 8^, wherein 8^ varies over the range [0, i;]. 

[0020] Fig. 3 illustrates a circular sensor array 30 comprising N sensors 8 distributed uniformly on a ring with radius 
R The N sensors are receiving a plane wave signal 22 arriving at an azimuth angle 9 and an elevation angle ^ as 
shown in Fig. 2. The A/ sensors 8 produce sensor time series 

50 

^ni^iX^-^ ^;i = l,...,/) 

55 indicative of the received plane wave signal 22. 

[0021] The plane wave signal response of the circular sensor array 30 is then expressed by: 



30 



35 
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5 X(f) Is a vector of the Fourier transform of the sensor time series 

10 

for a frequency bin / of interest. 

is a steering vector having its n^^ phase term for the plane wave signal 22 expressed by 

20 dXf. ^, ^) = oxp[j2 ^sin ^ cos(^ - ^)/ c] 

with e„ = 2 Tc/i/ A/, /7= 0,1,..., A/ - 1. The angles 9^ and e„are Indicating beam steering direction and sensor location, 
respectively. 

25 

30 is a diagonal matrix with off diagonal elements being zero and diagonal terms being weights of a spatial window to 
reduce side lobe structure. The spatial window is, In general, not unlfonn and depends on the sensor location and the 
beam steering direction. The beam power pattern 

Is then obtained by 

40 

p(f, /?. 4) = B(f, ^, 4y B'{f, i). 

An angular response of the circular sensor array 30 covers the range [0,2 tc] and, therefore, has no ambiguity with 

45 respect to the azimuth angle 6. 

[0022] Referring to Fig. 4 a decomposition process of a cylindrical array beamformer into coherent sub-sets of linear 
and circular array beamformers according to the Invention is shown. A cylindrical sensor array 40 as shown in Fig. 4 
comprises k sensors 8, wherein x = with A/ being the number of circular sensor arrays 30 and M being the 
number of sensors 8 on each circular sensor array 30, that is the number of line sensor arrays 20. Angular response 

^0 of the cylindrical sensor array 40 to a steered direction at 

55 

as indicated In Fig. 4 Is expressed as follows: 
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10 



15 



20 



25 



30 



35 



w^fjj is the (r, my^term of matrix W(Q, ^) comprising weights of a three dimensional spatial window to suppress sidelobe 
structures. Xf.fjj(f) Is the {r,m)^ temri of matrix X{f) expressing the Fourier Transfomi of a sensor time series of the m^^ 
sensor on the circular array. 

is a steering matrix having its (r, m)"' phase term for a plane wave signal expressed by 

drM' (l)= exp(72 /(r ^cos fi + Rsin ^cos(<?- ^))/ c), 

with R being a radius of the circular sensor array 30, 82 being a distance between subsequent circular arrays in 2- 

direction and Q„ = 2nm/ M, m=0,^ M - 1 indicating a position of a sensor 8 on the circular array 30. Inventor found 

a re-arranged fomn of equation (3), which is expressed as follows: 

b(/, <?. ^) = i; wiXC/". <?. f^i z ^,.™(/ d • (4) 

r-O L«-0 J 

dr(/* ^. ^) = exp(/2 ^ ^ cos ^ / c) 



is the r*'* temi of steering vector 



which is the kernel of a line array beamfonner as discussed above and ^Is the term of a spatial window for line an^ay 
40 spatial shading. 

dS/^ <?. i) = exp(;2 ^sin (^cos(^- ff„)/c) 

45 

is the m"* term of steering vector 

50 

which is the kernel of a circular array beamfonner as discussed above and Is the term of a spatial window for 
circular array spatial shading. Equation (4) expresses a cylindrical array beamformer comprising a product of two 
summations, the first being a linear array beamformer and the second being a circular array beamformer. According 
55 to the invention a beamfonning process for a cylindrical array as expressed by equation (4) is decomposed into two 
steps as shown in Fig. 4. In a first step circular an-ay beamfomiing is perfomned for each of the A/ circular sensor arrays 
30, each array comprising /Wsensors 8, to produce /Vbeam time series outputs. In a second step line array beamfonning 
is performed on the A/ beam time series outputs of the first step producing one beam time series for the steered direction. 
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The same process is then repeated for other beam steering directions of interest. Alternatively, the process is reversed 
perfomning line array beamforming first 

[0023] The decomposition process according to the Invention based on equation (5) substantially facilitates cylindrical 

array beamforming. The number of mathematical operations and the amount of memory required to perform these 
5 operations are substantially reduced by expressing the cylindrical beamformer as a product of two sums, Instead of a 
double summation as expressed by equation (3). It, therefore, eliminates the need for very large memory and very fast 
CPU's for real time system applications. This allows application of advanced beamforming algorithms for multidimen- 
sional arrays using currently available computer technology such as a conventional computer workstation in the com- 
pact processing unit 12. 

10 [0024] As is evident, the circular and line array beamformers resulting from the decomposition process according to 
the invention may be executed in parallel allowing optimal execution of the beamfomning operations using computer 
with parallel processors substantially enhancing performance of the beamforming process. This allows real time ap- 
plications of the ultrasound imaging system 100 using current available computer technology within the compact 
processing unit 12. 

15 [0025] Furthermore, because equation (4) is directly derived from equation (3) without further assumptions, the de- 
composed cylindrical array beamforming process according to the invention does not include a loss in the number of 
degrees of freedom. Therefore, the decomposed cylindrical array beamforming process according to the invention 
converges to an optimum solution as expressed by equation (3). 

[0026] Another advantage of the decomposition process according to the invention is a much simpler design and 
20 incorporation of three-dimensional spatial windows. A non-uniform shading window may be applied for each circular 
array beamformer to improve angular response with respect to azimuth angle G. A unifonn shading window may then 
be applied to the line array beamformer to improve angular response with respect to elevation angle <(>. Of course, other 
types of shading windows or combinations may be applied. 

[0027] Decomposition processes for planar and spherical arrays are very similar to the decomposition process for 
25 the cylindrical array described above. Detailed descriptions of these processes are published the inventor in A. Tawfik 
and S. Stergiopoulos: "A Generic Processing Structure Decomposing the Beamforming Process of 2-D & 3-D Arrays 
of Sensors into Sub-Sets of Coherent Processes", submitted to IEEE, J. Oceanic Eng., July, 1997. Of course, the 
decomposition according to the invention may also be applied to sensor arrays otherthan cylindrical, planar and spher- 
ical arrays. 

30 [0028] The beamforming process expressed by equations (1 ) and (2) for line sensor arrays 20 and circular sensor 
arrays 30, respectively, is a time delay beamfonning estimator being basically a spatial filter. However, optimum beam- 
forming requires beamforming filter coefficients to be chosen based on characteristics of noise received by the sensor 
array in order to optimize sensor array response. Algorithms for optimum beamforming using characteristics of noise 
received by the sensor array are called adaptive beamformers. Beamfomning filter coefficients of these algorithms are 

35 chosen based on a covariance matrix of correlated noise received by the sensor array. However, if the knowledge of 
the noise's characteristic is Inaccurate, perfomnance of the adaptive beamfonner will degrade significantly and may 
even result In cancellation of a desired signal. Therefore, it is very difficult to implement useful adaptive beamfomners 
in real time operational systems. Furthermore, for post processing such as matched filter processing the adaptive 
beamfonner has to provide coherent beam time series. In particular, the matched filter processing requires near-in- 

40 stantaneous convergence of the beamformer producing a continuous beam time series correlating with a reference 
signal. 

[0029] In adaptive beamforming, beamfomner response is optimized to contain minimum contributions due to noise 
and signals arriving from directions otherthan a direction of a desired signal. For the optimization, it Is desired to find 
a linear filter vector 

45 

so which is a solution to a constraint minimization problem that allows signals from a desired direction to pass with a 
specified gain. The minimization problem is expressed as follows: 

Minimize = Z ' (/, , ^R{f, )z(/, , ^ subject to T (/, , , = 1 . (5) 
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is the conventional steering vector based on equation (1). 



is a spatial correlation matrix of received sensor time series with elements 



wherein E{...) denotes an expectation operator and is sensor spacing between n^^ and sensor. The solution of 
the minimization problem is expressed by 



Equation (6) provides adaptive steering vectors for beamforming signals received by a W- sensor array. In frequency 
domain an adaptive beam at a steering angle 6^ is then defined by 



corresponding to conventional beams expressed by equation (1). 

[0030] In "Bearing Estimation of Multiple Broadband Sources using Steered Covariance Matrices", IEEE Trans. 
Acoust. Speech, Signal Proc, ASSP-37, pp. 1481-1494, 1989, J. Krolik and D.N. Swingler have shown that conver- 
gence time for a broadband source location Is reduced by using a space-time statistic called steered covariance matrix 
(STCM). This method achieves significantly shorter convergence times than adaptive algorithms based on a narrow- 
band cross spectral density matrix (CSDM) as taught by B. Van Veen and K, Buckley In: "Beamfomiing: a Versatile 
Approachto Spatial Filtering", IEEE ASSP Mag., pp. 4-24, 1988, without sacrificing spatial resolution. The number of 
statistical degrees of freedom available to estimate the STCM is approximately the product of time and bandwidth T 
X en/compared to an observation time (7"= M/ Fg, being a sampling frequency) in CSDM methods. This provides 
an improvement of approximately BW, being the size of a broadband source bandwidth, in convergence time. 
[0031] Output of a conventional beamfomier in frequency domain is expressed by equation (1). A corresponding 
output 




(6) 




(7) 




in time domain is then expressed as the weighted sum of steered sensor outputs 




(8) 



Since 
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IS an inverse fast Fourier transformation (IFFT) of 

5 



10 continuous beam time sequences are obtained from the output of a frequency domain beamformer using fast Fourier 
transformation (FFT) and fast convolution procedures. 

[0032] Using the beamformer output an expected broadband beam power 6(6) is given by: 

K^)=4^'n^1i} = h'E{r{t,,rS^X'n^M}f^' (9) 
wherein the vector h includes weights for spatial shading. The term 

20 

^t,,^)=E{x(t„ r„(^))} (10) 

is defined as the STCM In time domain and is assumed to be Independent of f/in stationary conditions. Supposing that 

25 



30 is the Fourier transfonn of the sensor time series and assuming that the sensor time series are approximately band 
limited a vector of steered sensor outputs 



is expressed by 



45 



is a diagonal steering matrix with elements Identical to the elements of the conventional steering vector 

50 



55 and is expressed as follows: 
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0 d,{f,,^ 



(12) 



10 



[0033] The STCM follows then directly from the above equations as 



15 



(13) 



20 



wherein the index k = IJ + 1 ..... / + H refers to frequency bins in a band of interest A/ and 



is the CDSM for the frequency bin 
25 [0034] In steered minimum variance algorithms (STMV) a spectral estimate of broadband spatial power is given by: 



30 



5(<?)=[7'cd(a/, <?)"'/]" 



(14) 



[0035] However, estimates of B(0) according to equation (14) do not provide coherent beam time series because 
they represent the broadband beam power output of an adaptive process. Therefore, according to the invention the 
estimation process of the STMV has been modified to detemfiine complex coefficients of 



35 



for all frequency bins in a frequency band of interest. 
40 [0036] The STIVIV algorithm is used in its original form to generate an estimate of <i>(Af, 6) for all frequency bands 
A/ of a received signal. Assuming stationarity across the frequency bins of a band A/ the estimate of the STMV Is 
considered to be approximately the same as a narrowband estimate 



for a center frequency /q of the band A/. Narrowband adaptive coefficients are then derived from 

50 
55 

Phase variations of 
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5 across the frequency bins / = /,/ + 1,..J + H are modeled by 

w,(/;.<^=exp[2/,T(A/,4], i = l,l + U.-,l + H, (16) 

10 

wherein ^^(Af, 6) Is a time delay term derived from 

«P,(Ar,^ = F[w.(A/.42^J. (17) 

Using the adaptive steering weights w„(^f, 9), provided by equation (16), adaptive beams are formed by 



B(f„^)=w(f„^)x(/,). (18) 

[0037] Fig. 5 shows schematically the processing steps of the STMV beamformer as discussed above. First the 
25 sensor time series are segmented, overlapped and transfonned into frequency domain using FFT. In a following step 
the STCM is detennined based on equations (10) and (13), The STCIVI is then inverted using Cholesky factorization 
and in order to determine an estimation of the adaptive steering vectors according to equation (15). The adaptive 
steering vectors are then Integrated In the frequency domain forming adaptive beams. Finally, adaptive beams In time 
domain are then formed through I FFT 

30 

bit,, <?. 4)= ifft{b(/„ 4' (19) 

35 

wherein overlap and concatenation of segments are discarded to form a continuous beam time series. 
[0038] Matrix inversion is a major issue for implementing adaptive beamfonners In real time applications. Standard 
numerical methods for solving systems of linear equations can be applied to solve for the adaptive weights. The nu- 
merical methods include: 

40 

• Cholesky factorization of the covariance matrix 




This allows the linear system to be solved by backsubstitution In terms of the received data vector. 

• QR decomposition of the received vector 

50 

m 

which includes conversion of a matrix to upper triangular form via rotations. The QR decomposition method has 
55 better stability than the Cholesky factorization but requires twice as much computational effort. 

• SVD (Singular Value Decomposition) method. The SVD method Is the most stable factorization technique but 
requires three time more computatlonai effort than the QR decomposition method. 
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[0039] For investigative studies of the beamforming process according to the present invention Cholesky factorization 
and QR decomposition techniques have been applied. No noticeable differences in perfornnance concerning stability 
have been found between these methods. Of course, for real time applications the fastest algorithm is pretended. 
[0040] Another major concern for implementing adaptive beamforming In real time applications Is the requirement 
of knowledge of second order statistics for the noise field, derived from equations (10) and (13). These statistics are 
usually not known but can be estimated from the received data by averaging a large number of Independent samples 
of the covariance matrix 

With /< being the effective number of statistically independent samples of 

and A/ being the number of sensors 8 the variance of the second order statistics for the noise field - the adaptive beam 
output power estimator detection statistic - Is Inversely proportional to f/C- A/+ 1). Theoretical estimations as well as 
empirical observations suggest that K has to be four times greater than N in order to get coherent beam time series 
at the output of the adaptive beamfonner. For arrays with a large number of sensors Implementation of adaptive beam- 
formers as statistically optimum beamformers requires averaging a very large number of Independent samples of 

which is not applicable in real time systems. A solution to this problem is updating of the covariance matrix using an 
exponential window in the form of a time varying adaptive estimation of 

at time t^^. The time varying adaptive estimation Is the exponentially time averaged estimator - geometric forgetting 
algorithm - at time t^^ : 



R" (fi) = (/,)+ 0 - /^^f.yifil (20) 

wherein |i is a smoothing factor (0<\l< 1) implementing exponentially weighted time averaging, 
[0041] Implementation of an adaptive beamformer with a large number of adaptive weights for a large number of 
sensors requires very long convergence periods eliminating dynamical characteristics of the adaptive beamformer to 
detect time varying characteristics of a received signal of interest. This limitation can only be avoided by reducing the 
number of adaptive weights. A reduction of the number of adaptive weights is achieved by introducing a sub-aperture 
processing scheme according to the present invention. 

[0042] Fig. 6 illustrates a sub-aperture configuration for a line array 20 of sensors 8. The line array is divided into a 
plurality of overlapping sub-arrays. In a first stage the sub-an^ays are beamfomied using a conventional beamformer 
generating a number of sets of beams equal to the number of sub-arrays for each steering direction. In a second stage 
adaptive beamfonning is perfomned on each set of beams steered In a same direction In space but each beam belonging 
to a different sub-array. A set of beams is equivalent to a line array consisting of directional sensors steered at a same 
direction with sensor spacing being equal to space separation between two contiguous sub-arrays and with the number 
of sensors being equal to the number of sub-arrays. 

[0043] Fig. 7 illustrates asub-apertureconflgurationforacircular array 30 with Msensors 8 according to the invention. 

A first circular sub-aperture consists of first M- G + ^ sensors with n = 1 ,2 M - G + 1 , wherein n is a sensor Index 

and G is the number of sub-apertures. A second circular sub- aperture consists of /If - G + 1 sensors with n = 2,3,..., 
/W - G + 2 . Sub-aperture formation Is continued until a last sub- aperture then consists oi M- G+^ sensors with n = 
G,G -f 1 M. In a first stage each circular sub-aperture is beamformed using a conventional beamformer generating 
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G sets of beams for each steering direction. The second stage includes adaptive beamforming on a set of beams 
steered in a same direction in space but each belonging to a different circular sub-aperture. Fig. 7 illustrates these 
sub-aperture steps for G = 3. 

[0044] Referring to Fig. 8, a sub-aperture configuration for a cylindrical array 40 according to the Invention is shown. 
The cylindrical array 40 comprises x = A//W sensors 8, wherein N \s the number of circular rings and M \s the number 
of sensors on each ring. Formation of sub-apertures is performed as follows: 

[0045] A first sub-aperture consists of the first ^A/ - G + 1) rings with n = 1 ,2 N - G+^ being a ring index and G 

being the number of sub-apertures, in each ring a first set of /If - G + 1 sensors, with m = 1 ,2,... M - G + ^ being a 
sensor index for each ring, is selected. The first sub-aperture forms a cylindrical array cell 42, Indicated by small circles 
in the upper right picture of Fig. 8. Conventional beamforming is then applied to the cylindrical array cell 42 using the 
decomposition process according to the invention as described above. For a given steering direction comprising azi- 
muth angle Og and elevation angle ^gihe multidimensional sub-aperture beamformer provides beam time series 

with subscript g= 1 Indicating the first sub-aperture. 

[0046] A second sub-aperture consists of a subsequent set of (A/ - G+ 1) rings with n = 2,3 N- G+ 2. In each ring 

a consecutive set of - G + 1 sensors with m = 2,3,...M - G + 2 is selected. These selected sensors then fonri a second 
sub-aperture cylindrical array cell. Again conventional beamfonning is then applied and for a given steering direction 
the multidimensional sub-aperture beamformer provides beam time series 

with subscript g=2 indicating the second sub-aperture. 

[0047] Sub-aperture fomriation is continued until a last sub- aperture then consists of a subsequent set of (A/ - G + 

1) rings with n = G, G + 1 N . In each ring a last set of Af - G + 1 sensors with m = G, G + 1,..., M is selected. 

Conventional beamforming is then applied and for a given steering direction the multidimensional sub-aperture beam- 
former provides beam time series 

with subscript g= G indicating the last sub-aperture. 

[0048] A second stage of beamforming comprises an adaptive beamfomier on a line array consisting of, for example, 

G = 3 beam time series 

9 = 1 ,2 G, as shown in Fig. 8. For a given pair of azimuth and elevation steering angles 

the cylindrical adaptive beamforming process is reduced to an adaptive line array beamformer. The adaptive line an^ay 
beamfonner comprises only three beam time series 
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g = 1 ,2,3 with spacing 

between two contiguous sub-aperture cylindrical cells, wherein (R2 n/ M) \s the sensor spacing on each ring and 
wherein 5 is the distance between each ring along the z-axis of the cylindrical array 40. The adaptive line array beam- 
former provides one or more adaptive beam time series with steering centered on the pair of azimuth and elevation 

steering angles 

[0049] Because of the very small number of degrees of freedom in each sub-aperture the adaptation process expe- 
riences near-instantaneous convergence. Furthermore, the multidimensional sub-aperture beamforming process ac- 
cording to the invention may include a wide variety of adaptive noise cancellation techniques such as MVDR and GSC 
as discussed in: 

S. Stergiopoulos: "Implementation of Adaptive and Synthetic Aperture Processing Schemes in Integrated Active- 
Passive Sonar Systems", Proceedings of the IEEE, 86(2), pp. 358-397, Febmary, 1998. 
B. Windrow et al.: "Adaptive Antenna Systems", Proceedings IEEE, 55(12), pp. 2143-2159, 1967; 
N.L. Owsley: "Sonar Array Processing", 8. HaykIn, Editor, Prentice-Hall Signal Processing Series, A.V. Oppenhelm 
Series Editor, pp. 123, 1985; and, 

B. Van Veen and K. Buckley: "Beamfomriing: a Versatile Approach to Spatial Filtering", IEEE ASSP Mag., pp. 4-24, 
1988. 

[0050] Furthermore, the sub-aperture configuration according to the invention is applicable to other multidimensional 
arrays such as planar arrays and spherical arrays. Decomposition, sub-aperture formation as well as implementation 
of adaptive beamfomners for planar and spherical arrays are similar to corresponding steps for cylindrical an-ays dis- 
cussed above. In particular, for planar arrays, formation of sub-apertures is based on the sub-aperture adaptive concept 
for tine arrays as illustrated in Fig. 9. Similarly, formation of sub-apertures for spherical arrays is based on the sub- 
aperture adaptive concept of circular arrays. 

[0051 ] Fig. 1 0 Illustrates schematically the modular stnjcture of the signal processing flow for a 3D cylindrical array 
of an ultrasound imaging system according to the invention. Line and circular array beamfomners are implemented as 
Finite Impulse Response (FIR) filters. Reconfiguration of the different processing modules in Fig. 1 0 allows application 
of the signal processing flow to a variety of ultrasound imaging systems with, for example, line, planar or spherical 
sensor arrays. 

[0052] A first processing module 50 includes partitioning of time series received from a sensor array, computation 
of an initial spectral FFT of the partitioned time series, selection of a frequency band of interest via band-pass FIR 
filters and downsampling. The output of module 50 comprises a continuous time series at a reduced sampling rate. 
[0053] Second module 52 and third module 53 comprising circular array and line array beamformers provide con- 
tinuous directional beam time series using a FIR scheme for spatial filtering via circular convolution. Segmentation and 
overlap of an input of the beamformers takes care of wraparound errors due to fast-convolution signal processing 
operations. The overlap is equal the effective FIR filter's length. 

[0054] IVIoduie 54 is used for processing echoes obtained from active ultrasound applications. 
[0055] IVIoduie 56 includes final processing steps of a temporal spectral analysis. 

[0056] Finally, data normalization processing schemes are used to map output results Into a dynamic range of display 

devices. 

[0057] In the module 56 use of verniers and the temporal spectral analysis - incorporating segment overlap, window- 
ing and FFT coherent processing - provide narrow band results for all the beam time series. Normalization and OR- 
ing are final processing steps before displaying the output results. 

[0058] In the active module 54 application of a matched filter - or replica correlator - on the beam time series provides 
broadband processing. This allows detection of echoes as a function of range and bearing for reference waveforms 
transmitted by active transducers of ultrasound systems. 

[0059] Referring to Fig. 11. a schematic diagram of a signal processing flow for the sub-aperture configuration of 
multidimensional arrays according to the invention is shown. Module 60 includes fomiation of sub-apertures according 
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to the invention. In a next step, sensor time series of each sub-aperture are beamformed in module 62 using a con- 
ventional beamformer for multidimensional arrays according to the invention as shown In Fig. 10 providing a beam 
time series for each sub-aperture. A second stage of beamforming Is performed in module 64. The second stage 
Includes adaptive beamforming of the beam time series obtained in module 62. The adaptive beamfonner provides 
5 one or more adaptive beam time series with steering centered on same azimuth and elevation angles as the conven- 
tional beams. 

[0060] The signal processing flow for the sub-aperture configuration of multidimensional arrays as shown in Fig. 11 
may be incorporated in modules 52 and 53 of the signal processing flow shown in Fig. 10. 

[0061 ] In another embodiment of signal processing according to the invention the sub-aperture adaptive beamformer 
10 is combined with a synthetic aperture beamformer - ETAM algorithm - for further improvement of the image resolution. 
The synthetic aperture beamfonner provides sensor time series corresponding to a synthesized sensor array compris- 
ing a larger number of sensors than the real • physical - sensor array used for providing input sensor time series to the 
synthetic aperture beamfonner. Details concerning the synthetic aperture beamformer are disclosed in US Patent No. 
4,930,111 issued to Sullivan E.D. and Sterglopoulos S. in 1990. In a first step the sensor time series are processed 
15 using the synthetic aperture beamfonner and are provided as input to the sub-aperture adaptive beamformer. Alter- 
natively, this method may be used to obtain a satisfying image resolution when smaller sensor arrays are deployed in 
compact ultrasound imaging systems. 

[0062] The sub-aperture adaptive beamformer according to the Invention has been tested with simulated as well as 

with real data sets In order to demonstrate that it achieves near instantaneous convergence and, therefore, provides 
20 output beamtime series having sufficient temporal coherence and correlate with the reference signal. Provision of 
output beamtime series having sufficient temporal coherence is essential for matched filter processing or temporal 
spectral analysis in passive and active ultrasound imaging systems. Referring to Fig. 1 2, the power of beam response 
to active CW pulses as function of steering angle for various types of beamfomiers - conventional beamformer on a 
48 sensor planar array (1); sub-aperture adaptive beamformer according to the Inventionon a 48 sensor planar array 
25 (2); synthetic aperture (ETAM algorithm) beamfonner extending a physical 48 sensor planar array into a synthetic 1 44 
sensor planar array (3); combined synthetic aperture and sub-aperture adaptive beamformer according to the invention 
on a 48 sensor planar array (4); and. conventional beamformer on a 144 sensor planar array (5) - is shown. As is 
apparent from the results presented in Fig. 12, the sub-aperture adaptive beamformer according to the invention pro- 
vides an angular resolution similar to a conventional beamformer applied on a three times larger planar array. Further 
30 improvements to angular resolution are obtained by combining the sub-aperture adaptive beamformer with the synthetic 
aperture beamfonner according to the invention as shown in Fig. 1 2 at the expense of more computational effort being 
involved in the beamforming process, 

[0063] Furthennore, the sub-aperture adaptive beamformer according to the Invention has been tested with synthetic 
data sets including broadband FM pulses with a repetition rate of a few minutes to demonstrate the efficiency of the 

35 sub-aperture adaptive beamformer according to the invention to achieve near-instantaneous convergence In order to 
provide coherent beamtime series for the broadband FM pulses. This Is necessary to test, because adaptive processing 
schemes require at least a few iterations to converge to an optimum solution. The output beamtime series have been 
processed using a matched filter and their arrangement as GRAMs provides a waterfall display of ranges (depth) as 
a function of beam-steering as shown In Fig. 13. For each beamfonner the horizontal axis includes 20 steering beams 

40 In a 10-degrees angular sector centered in a -20 degree lool<-direction. The vertical axis represents time delay or depth 
penetration of the signal within a human body. Thus, the detected echoes along the vertical axis of each window shown 
In Fig. 13 represent reflections from simulated objects such as organs. From left In Fig. 13 the different beamformers 
used are; a conventional beamfonner on a 48 sensor array; a sub-aperture adaptive beamformer according to the 
inventionon on a 48 sensor array; a synthetic aperture (ETAM algorithm) beamformer extending a physical 48 sensor 

45 array into a synthetic 1 44 sensor array; a combined synthetic aperture and sub-aperture adaptive beamformer accord- 
ing to the invention on a 48 sensor array; and, a conventional beamformer on a 144 sensor array. As is obvious, the 
conventional beamformer applied on the 48 sensor array - first window on the left in Fig, 13 - has failed to resolve the 
two echoes. This represents the performance characteristic of current state-of-the-art ultrasound imaging systems. 
The following three windows in Fig. 13 show results of the adaptive beamformers according to the invention applied 

50 to a 48 sensor array. As is evident, in all three windows two echoes are clearly visible. The resolution of the adaptive 
beamformers according to the Invention Is comparable to the resolution of a conventional beamformer applied on a 
three times larger sensor array as shown in the window on the right in Fig. 13. 

[0064] In order to demonstrate the improvement, the adaptive beamformer according to the present invention has 

been implemented within a general purpose ultrasound imaging system for medical diagnostic applications as shown 
55 In Fig. 1 4. The system comprises: a line array 20 of 32 sensors 8; a state-of-the-art signal processing unit 72; a custom- 
made A/DC 74 with 12 bit and 45MHz sampling frequency per channel; and, a computer workstation 76 for adaptively 
beamfomn the sensor time series provided by the line array 20 and for image reconstruction and for 3D visualization. 
[0065] Fig. 15 shows typical 3D Images of a fetus' skull. The two top Images and the bottom left image have been 
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obtained using state-of-the-art ultrasound systems and signal processing whereas the bottonn right image has been 
obtained by applying the adaptive beamformer according to the present invention to the signals captured by a state- 
of-the-art ultrasound system. As Is evident, application of the adaptive beamformer according to the present Invention 
substantially enhances the image resolution. 
5 [0066] Further enhancement of the Image resolution is obtained by applying the adaptive beamformer according to 
the present invention to multidimensional - 2D or 3D - sensor arrays allowing better angular resolution of the signals 
during the beamforming process. 

[0067] As shown by the results presented In Figs. 1 2, 1 3 and 1 5 the adaptive beamfomier according to the present 

invention achieves near-instantaneous convergence and enhances image resolution substantially. Therefore, the adap- 
10 tive beamformer is highly advantageous for applications in ultrasound imaging and may be implemented into existing 
ultrasound systems for Improving image resolution. In a preferred embodiment the adaptive beamformer according to 
the invention is used to provide a compact field-deployable ultrasound system for providing high resolution images in 
real time. This is achieved by applying the adaptive beamfomier on signals obtained from compact multidimensional 
arrays and using a sub-aperture structure reducing the degrees of freedom within each sub-aperture thus minimizing 
15 the computational effort for beamforming multidimensional arrays. Therefore, the beamforming process requires only 
a conventional workstation and is perfonned in real time. 

[0068] Of course, numerous other embodiments may be envisaged without departing from the spirit and scope of 
the invention. 

20 

Claims 

1. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system comprising the steps of: 

decomposing the multidimensional beamformer into two coherent subsets of circular array beamfomners and/ 
or line an-ay beamformers, a first subset comprising subsequent beamformers in a first coordinate direction 
of the multidimensional array and a second subset comprising subsequent beamformers in a second other 
coordinate direction of the multidimensional array; 

beamforming for a predetermined beam steering direction of data relating to the sensortime series by applying 
the subsequent beamformers of the first subset, each beamfomier producing a beam time series; and, 
beamforming for the predetermined beam steering direction each beam time series of the first subset of beam- 
fonners applying the subsequent beamfomners of the second subset for the steered direction producing one 
beam time series for the beam steering direction. 

2. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1 , wherein the multidimensional sensor array is 
a planar array. 

40 3. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1 , wherein the multidimensional sensor array is 
a cylindrical array. 

4. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
45 sional array of an ultrasound imaging system as defined in claim 1 , wherein the multidimensional sensor an^ay is 

a spherical array. 

5. A method for multidimensional beamfonning sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1 , wherein the subsets of beamformers of each 

50 coordinate direction are executed in parallel. 

6. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1 , comprising weights of a three dimensional 
spatial window to suppress sidelobe structures. 

55 

7. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1, wherein the multidimensional beamformer 
comprising a double summation of terms in two different coordinate directions of the multidimensional anray Is 
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transformed into a product of two summations, each summation comprising terms in one coordinate direction of 
the multidimensional array. 

8. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1 , wherein the beamformers of the subsets are 
adaptive beamfomiers. 

9. A method for multidimensional beamfonning sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1 , wherein non-unifomri shading windows are 
applied to the beamformers of a subset of beamfomners. 

10. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1 , wherein the line and circular array beamformers 
are implemented as Finite Impulse Response (FIR) filters for spatial filtering via circular convolution. 

11. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system as defined in claim 1 , comprising broadband processing of the beam 
time series for each steering direction using a matched filter for detection of echoes as a function of range and 
bearing. 

12. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system using a sub-aperture configuration, the method comprising the steps 
of: 

a) dividing the multidimensional beamfonner Into a plurality of subsequent sub-apertures; 

b) decomposing each sub-aperture into two coherent subsets of circular array beamformers and/or line array 
beamformers, a first subset comprising subsequent beamformers in a first coordinate direction of the multidi- 
mensional array and a second subset comprising subsequent beamfomners in a second other coordinate 
direction of the sub-aperture; 

c) conventional beamfonning each sub-aperture for a predetemnined beam steering direction a Fourier trans- 
form of the sensor time series by applying the subsequent beamformers of the first subset each beamfonner 
producing a beam time series; 

d) conventional beamfonning each sub-aperture for the predetermined beam steering direction the beam time 
series produced by the beamfomners of step c) by applying the subsequent beamformers of the second subset 
forthe steered direction producing one beam time series for the beam steering direction for each sub-aperture; 

and, 

e) adaptive beamforming on line arrays, each line array comprising beam time series of different subsequent 
sub-apertures In one coordinate direction, providing one or more beam time series for the beam steering 
direction. 

13. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system using a sub-aperture configuration as defined in claim 12, wherein 
the multidimensional sensor array is a planar array 

14. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system using a sub-aperture configuration as defined in claim 12, wherein 
the multidimensional sensor array Is a cylindrical array 

15. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system using a sub-aperture configuration as defined in claim 12, wherein 
the multidimensional sensor array is a spherical array. 

16. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system using a sub-aperture configuration as defined In claim 1 2, the method 
comprising MVDR adaptive noise cancellation, 

17. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system using a sub-aperture configuration as defined in claim 1 2, the method 
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comprising GSC adaptive noise cancellation. 

18. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system using a sub-aperture configuration as defined in claim 12, wherein 

5 the adaptive beamfonner is combined with a synthetic aperture beamfomier 

19. A method for multidimensional beamforming sensor time series provided by sensors deployed in a multidimen- 
sional array of an ultrasound imaging system using a sub-aperture configuration as defined In claim 12, wherein 
the snesor time series are beamformed in real time. 

10 

20. A method for beamforming sensor time series provided by sensors of an ultrasound imaging system using a co- 
herent broad band adaptive beamformer, the method comprising the steps of: 

segmenting the continuous sensor time series into a set of overlapped data sets; 
15 calculating a FFT of each overlapped data set producing a set of Fourier transforms of the overlapped data 

sets for different frequency bins; 

forming a cross spectral density matrix from the Fourier transfonns of the overlapped data sets for each fre- 
quency bin and each predetermined steering direction; 

fomiing a steering covariance matrix using the cross spectral density matrix and a diagonal matrix of conven- 
20 tional steering vectors, one steering covariance matrix for each steering direction and a frequency band of 

interest; 

inverting the steering covariance matrices; 

estimating adaptive steering vectors by assuming stationarlty across frequency bins of a frequency band of 
interest and considering an estimate of the steering covariance matrix being the same as a narrow band 
25 estimate for a center frequency of the frequency band of interest; 

determining narrow band adaptive steering weights using the estimate of the adaptive steering vectors; 
fomning adaptive beams in frequency domain from the Fourier transfomn of the overlapped data sets and the 
adaptive steering weights; 

fonning adaptive beams in time domain through I FFT; and, 
30 determining continuous beam time series by discarding overiap and concatenation of segments. 

21. A method for beamforming sensor time series provided by sensors of an ultrasound imaging system using a co- 
herent broad band adaptive beamfomier as defined in claim 20, wherein the steering covariance matrices are 
inverted using Cholesky factorization. 

35 

22. A method for beamfonning sensor time series provided by sensors of an ultrasound imaging system using a co- 
herent broad band adaptive beamformer as defined in claim 20, wherein the steering covariance matrices are 
inverted using QR decomposition. 

40 23. A method for beamforming sensor time series provided by sensors of an ultrasound imaging system using a co- 
herent broad band adaptive beamformer as defined in claim 20, wherein the steering covariance matrices are 
inverted using SVD (Singular Value Decomposition) method. 

24. A method for beamfomning sensor time series provided by sensors of an ultrasound imaging system using a co- 
45 herent broad band adaptive beamformer as defined In claim 20, wherein the steering covariance matrices are 

updated using a time varying adaptive estimation, the time varying adaptive estimation being an exponentially time 
averaged estimator. 

25. A field deployable 3D ultrasound imaging system for producing high resolution 3D images of an object In real time, 
50 the ultrasound system comprising: 

a source for emitting ultrasound waves; 

a compact adaptive multidimensional sensor array for capturing reflections of the ultrasound waves, the ul- 
trasound waves being reflected by different structures within the object, and for providing sensor time series 
55 indicative of the reflected ultrasound waves; 

a compact processing unit for: 

receiving the sensor time series produced by the multidimensional sensor an^ay; 
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processing the sensor time series In order to produce continuous beanntime series by: 

decomposing a multidimensional beamformer Into sub-apertures comprising coherent subsets of cir- 
cular array beamformers and/or line array beamfonners; 
conventional beamfomiing circular arrays; 
adaptive beamfonning line arrays; 



reconstructing 3D Images from the beam time series In real time; 

and, 

a display for displaying the reconstructed 3D images in real time. 

26. A field deployable 3D ultrasound imaging system for producing high resolution 3D images of an object in real time 
as defined in claim 25, wherein the compact adaptive multidimensional sensor array comprises a planar array 

27. A field deployable 3D ultrasound imaging system for producing high resolution 3D images of an object in real time 
as defined In claim 25, wherein the compact adaptive multidimensional sensor array comprises a cylindrical array. 

28. A field deployable 3D ultrasound imaging system for producing high resolution 3D images of an object in real time 
as defined In claim 25, wherein the compact adaptive multidimensional sensor array comprises a spherical array. 

29. A field deployable 3D ultrasound imaging system for producing high resolution 3D images of an object in real time 
as defined In claim 25, wherein the compact processing unit comprises a conventional computer workstation. 

30. A field deployable 3D ultrasound imaging system for producing high resolution 3D images of an object in real time 
as defined In claim 28, wherein the conventional computer workstation comprises a plurality of processors for 
parallel processing beamfonners of a subset. 
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Fig.l 
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Figure 6 
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Figure 8 
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APAPTIYE BEAMFORMER 
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Figure 11 
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Figure 13 
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